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PREFACE  AND  SUMMARY 

I 

One  aspect  of  RAND  work  on  the  development  of  com¬ 
puter  models  to  represent  physiological  subsystems  is  con¬ 
cerned  with  the  resolution  of  attendant  problems  of  theo¬ 
retical  and  computational  mathematics.  In  a  previous 

t 

Memor andum- -Mass  Action  Laws  and  the  Gibbs  Free  Energy 
Function >  N.  Z.  Shapiro  and  L.  S.  Shapley,  RM-3935-PR, 

I 

December  196!}^-- the  use  of  the  Gibbs  free  energy  function 
vis-a-vis  the  mass  action  laws  for  the  analysis  of  complex 
chemical  systems  was  considered.  The  present  Memc randum 
continues  that  discussion,  obtaining  an  inequality  concern 
ing  the  effect  on  a  chemical  system  of  changing  its  free 


energy  parameters. 
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INTRODUCTION 


In  this  Memorandum  we  obtain  an  inequai ity  which 
furnishes  information  about  the  behavior  of  a  chemical 
equilibrium  system  when  its  free  energy  parameters  are 
changed. 

The  result  is  proved  by  viewing  the  solution  of  a 
chemical  equilibrium  problem  as  the  result  of  minimizing 
the  free  energy  function.  Were  the  solution  viewed  in 
terms  of  the  mass  action  laws,  the  proof  would  become 


much  more  difficult. 
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ON  THE  BEHAVIOR  OF  A  CHEMICAL  EQUILIBRIUM  SYSTEM 
WHEN  ITS  FREE  ENERGY  PARAMETERS  ARE  CHANGED 

Theorem  1.  Let  P  and  be  chemical  equilibrium 
•k 

problems  which  are  identical  except  possibly  for  their 
free  energy  parameter  vectors,  c  and  c^,  (That  is,  P  and 
P^*^  contain  the  same  species,  compartmentalized  in  the 
same  way  and  obeying  identical  mass  balance  laws.  Ax  =  b* 
In  other  words,  is  the  problem  obtained  from  P  by  re¬ 
placing  c  by  ,)  Let  x  be  any  solution  of  P  and  let 
be  any  solution  of  P^^.  (Note:  We  do  not  assume  that  P 
or  P^^  necessarily  has  a  unique  solution.  Thus,  x  may  be 
one  of  possibly  many  solutions  of  P.)  (Note:  We  do  not 
assume  that  x  or  is  positive;  that  is,  some  of  their 
components  may  be  zero.)  It  then  follows  that 


X 


•  (c 


X-  (c 


(1) 


Furthermore,  if  c  -c  is  not  representable  as  a  linear 
combination  of  the  rows  of  the  constraint  matrix  A,  and 


•k 

For  a  detailed  introduction  to  the  subject  matter  of 
this  Memorandum,  see  Part  I  of  RM-3935-PR,  Mass  Action  Laws 
and  the  Gibbs  Free  Energy  Function,  N.  Z.  Shapiro  and  L.  S. 
Shapley,  December  1963,  and  also  the  publications  referenced 
in  that  Memorandum. 
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if  all  the  components  of  either  x  or  of  x 
and  x^)  are  positive,  then 


(or  of  both  X 


X 


-c) 


-c)  . 


Before  proving  Theorem  1,  we  give  two  examples  of 
applications  of  the  theorem. 

Example  1.  If  a  chemical  equilibrium  problem  with  a 
solution  X  is  altered  by  increasing  one  of  the  free  energy 
parameters  c^  ,  and  if  x^^  is  any  solution  of  the  resulting 
problem,  then 


Furthermore,  if  x  is  positive  and  if  x^  is  not  de¬ 
termined  solely  by  the  mass  balance  constraints,  then 

#  , 

X  .  <  X  .  . 

J  J 

Example  2.  Let  a  chemical  equilibrium  problem,  having 
a  positive  solution,  x,  be  altered  by  increasing  each  of 
the  c.'s  in  one  of  the  compartments  by  the  same  positive 


I 
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imount.  (Note:  This  co-responds,  loosely  speaking,  to 

increasing  the  pressure  of  that  compartment.)  Let  be 

any  solution  of  the  resulting  problem.  Let  x  and  denote 

the  sum  of  the  x.'s  and  of  the  xf ' s  in  the  affected  com- 

J  J 

partment.  Then,  if  x  is  not  determined  solely  by  the  mass 
balance  constraints,  we  have 

- 

X  '  X  . 

Theorem  1  follows  from  a  more  general  result.  To 
state  this  more  general  proposition:  Let  n  be  a  positive 
integer.  Let  be  Euclidian  n-dimensiona  1  space.  Let 
H  be  a  convex  subject  of  E*\  Let  be  the  interior  of 
H  relative  to  the  linear  nuinifold  on  E^^  generated  by  H. 

Let  F  be  a  real-valued  convex  function  defined  on  H. 

If  u  is  any  element  of  H  and  if  9  is  any  element  of 
E^  such  that  u  +  t9  is  an  element  of  H  for  all  sufficiently 
small  positive  real  t,  then  the  directional  derivative 
F^(u)  is  defined  by 

F'(u)  =  lin  — - J - ^ 

^  t-(H-  ^ 


(2) 
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It  is  well  known  that  the  convexity  of  F  assures  the 
existence  (if  we  accept  -  as  a  limit)  of  Fq(u).  It  is 
also  well  known  that  if  u  and  v  are  elements  of  H,  then 

F(v)  -  F(u)  F^_^(u)  .  (3) 

We  will  say  that  F  is  weakly  differentiable  at  a 
point  u  in  H,  if  for  all  0  (for  which  u  +  t0  is  in  H  for 
all  real  t  with  sufficiently  small  absolute  values)  for 
which  FQ(a)  -  0  and  F^q(u)  ^  0,  we  have 


F^(u)  =  FIq(u)  =  0  .  (4) 

(Note;  Weak  differentiability  follows  from  the  stronger 
assumption  that  F^q(u)  =  -F^(u)  for  all  0  (for  which 
u  +  t0  is  in  H  for  all  t  having  sufficiently  small  ab¬ 
solute  values),) 

Let  A  be  any  element  of  E^.  Let  F^^  be  a  real-valued 
function  defined  on  H  by  means  of 


F^\u)  =  F(u)  +  4»u  , 


(5) 
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Then,  Theorem  1  will  follow  from: 

Theorem  2.  Let  x  be  a  member  of  the  minimum  set  of 
F  on  H,  and  let  be  a  member  of  the  minimum  set  of 
on  H,  rnen, 


A-  X 


(6) 


Furthermore,  if; 

(A)  F  is  weakly  differentiable  at  x 

(B)  X  is  in 

(C)  A«u  is  not  constant  for  all  u  in  H 


or  if: 


(/) 

a*) 

(C*) 


F  is  weakly  differentiable  at  x 

#  .  .  uO 

X  IS  in  H 

A»u  is  not  constant  for  all  u  in  H 


then; 


<  A-  X 


(7) 


Proof  of  Theorem  2.  Let  us  first  observe  that  for  any  u 
in  H  and  for  any  0  for  which  u  +  t0  is  in  H  for  all  suf¬ 
ficiently  small  positive  t 
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F*(u)  =  F^(u)  +  A-0  . 


(8) 


In  particular, 


(x)  =  (x)  +  A-(x^'^-x)  .  (9) 

X  -X  X  -X 


But,  since  x  is  in  the  minimum  set  of  F,  we  must  have 


F'y^  (x)  0  .  (10) 

X  -X 


But  by  (3) 


^  F^'^  (x)  i  (By  (9)  and  (10))  d-  (x''-x) 

X  -X 

(11) 


# 


# 


But,  since  x  is  in  the  minimum  set  of  F  ,  we  have 


#V 

F  (x  ) 


F*(x) 


(12) 


then,  (6)  follows  (11)  and  (12), 

To  prove  that  (7)  follows  from  (A),  (B)  ,  and  (C) ,  or 
from  (A^^)  ,  (B^^),  and  (C^^)  ,  it  would  be  sufficient  (by  the 


symmetry  obtained  when  A  is  replaced  by  -A)  to  show  that 


(7)  follows  from  (A),  (B)  ,  (C). 


Let  us,  therefore,  assume  (A),  (B) ,  and  (C) ,  and  let 
us  also  assume  that  (7)  is  false.  But,  since  we  already 
proved  (6),  we  have  that 


A-  (x^^-x)  =  0  . 


(13) 


But  this  and  (11),  then,  yield 

f"(x)  .  (14) 

so  that  (12)  and  (14)  yield 

it  ii 

F^(x  )  =  F  (x)  .  (15) 

Let  u  be  a  point  in  H  guaranteed  by  (C)  such  that 
A*u  y  ^‘X.  Let  0  =  u-x.  By  (B) ,  x  +  t0  is  in  H  for  all 
t  with  sufficiently  small  absolute  values. 

Since  x  is  in  the  minimum  set  of  F,  we  know  that 


Fq(x)  •  0,  and  F]^q(x)  0.  Hence,  by  (A) 
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Hence,  by  (8), 


Fq(x)  =  A*G  and  ^{x)  =  -AO 


Since  ^-0  ^  0,  we  have  (if  necessary  by  replacing  0  by  -0) 


# ' 

Fe(x)  -  0  . 


(16) 


But,  since  x^^  is  in  the  minimum  set  of  (15)  show  i 
that  X  is  in  the  minimum  set  of  ,  a  fact  which  contra¬ 
dicts  (16).  QED. 


